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INTRODUCTION 


Rpsearch ot Ohio University is currently being done to develop a low-cost navigation 
package utilizing the Loran-C hyperbolic lattice system. One of the components of this 
package is a receiver to pick up the Loran-C signal* und present them in the form of time 
differences. The time differences do not provide much information by themselves; it is 
necessary ta convert them to a more useful form, such as geographic coordinates or range- 
bearing Information. Although this has been done with charts and maps, today's high speed 
and inexpensive computing systems make the use of Loran-C much more attractive for aviation 
Installations. This report will provide some of the simpler mathematical equations which may 
be used In Loran-C navigation calculations. Specifically, a technique will be presented to 
allow Loran-C time differences to be predicted at a location. This is useful for receiver 
performance work, and as a tool for more complex calculations, such as position fixing. 


II REVIEW OF LORAN-C SYSTEM 

The standard Loran-C configuration consists of a master station and two or more 
secorsdary "slave" stations, synchronized to the master (see Figure I). At the start of the 
group repetition interval (GRI), the r ster station tranTniis o set of pulses which are 
received at the secondary stations after an approprio're time delay. Each secondary then 
inserts a unique coding delay before re-tronsmitting the pulses. Reception of o master- 
secondary pair then defines a line-of-position (LOP) which, if plotted for a constant time 
difference (TD) reading, results In a hyperbola with the master and secondary as foci 
Reception of two master-secondory pairs result in two TD readings corresponding to the 
intersection 'f two LOP's. Conversions rrxjy now be done to obtain a more meaningful 
position fix. Using only two LOP's also results in a second ambiguous position, although 
this ordinarily does not present a serious problem in position fixing. The reception of 
Loran-C stations and geii“ration of time difference information is o classic problem 
and will not be considered in this text. 


Ill CALCULATION OF TIME DIFFERENCES 

The receiver meosures the time difference between the arrival of the master signal 
and the arrival of the secondary signal. This is sometimes known os the difference mode. 
Since the receiver cannot directly measure the distance to each of the stations, it effec- 
tively measures the time difference, or length, from the master to the secondary, plus the 
secondary coding delay, plus the distance from secondary to receiver, minus the distance 
from receiver to master. In equation form: 

TD= B + At-T-T 
s m 

where B represents the baseline length from master to secondary, Ais thr coding delay 
at the secondary, and T^ and T^ are the one-way baseline lengths from the receiver to the 

secondary and master respectively. The fixed quantities p and Adepend on the system 
geometry and the two baseline times are the unknowns which must be calculated. 
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A geodesic model of the earth will now he developed to calculate distances between 
two points. The following discussion will he done in conjunction with Figure 2, A sphere 
with radius a Is circumscribed around the best-fitting ellipsoid of the earth. The angle d> 
(latitude) is found from a segment perpendicular to the tangent at a point P on the ellipsoid. 
Angle P is found from the segment joining the vertical projection of P on the spheroid. 

From Figure 2: 

X » a cos p 

y » b sin p 


and 


tan d> » -dv''dy 

substituting 

tan P » (b a) tan<t». 

If two points having coordinates (<b. , \ ) ui:d , \^) arc considered, the circular 
- distance is; ^ ^ 


arc 


cos X = sin pj sin P^ cos P^ cos cos (X^ - X^) 

By multiplying the arc distance just obtained by the semi-major ladius of the earth, the 
distance may be obtained in the same units used for the radius (miles, kilometers, micio- 
secorsds, etc.). Solving tl«e distance on the geold proceeds as follows 

P “ (a-b) (x - sin x) / 4(1 + cos x) 


Q = (a-b (x t- sin x) 4(1 - cos x) 
Aq = (sin Pj sin 

2 

Bp = (sin P^ - sin P^) 

Finally 


d - D - AqP - BgQ 

where d is the geodesic arc length and D is the spherical arc length (in the same units), 


By following the above procedure, the arc length may he found for receiver-to- 
master and receiver-to-secondaiy and used to solve the difference equation given earlier. 

Figure 3 shows the a!:ove procedure os implemented in a computer program using tlie stan- 
dard FORTRAN-IV progromming languoge. The subprogram "DIS" calculates the geodesic 
arc length from the coordinates of two points. Figure 4 gives an explanation of the variables 
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pile: time 


FOUTRftV c 


OHIO UMVEPSITY AVIONICS ENGINFfRING CENTER 


IMPLICIT PEAl.CL) 

OImcmsION CCRD(6)tSTATC5t?),9ETA(A)tOELTMA),TO<A),TC5) 
COMMON £>hIR,LAMR 


TIMOOOlO 

timooc?o 

TIM00050 


CATA STAT/0.59A50 74 3*C.<*71ttl301.0.8163e789,n,7200037l«0.6955A99«l.TIMCOO<»C 


16 

t 


\ 339fi 391 9, 1 ,39P?6820» 0.92006991,1 .221333921 f 1.52693116/ 
DATA RE T A /2b95.51 ♦0309.66,359 1.31 ,3560 .71/ 

DATA DELTA/1.1E9,2.0E9,9.9E9,6.5E9/ 

DATA ROS/1 .79532925E-2/ 

R0INT 6 

PORmak *OENTER LATITUDE AND LONGITUDE CP POSITION; ••CR»* 

:nate.»/» sddd '•m ss.sm 

PEA0(5,l,ENDsl5)C0>»'^ 


1 PO0MAT«P9.O,lif,P2.O,l<,f9.O) 

PhI'> = POS»(COPD(l)»<CORO<2)*CORDI 3)/60. )/60. ) 
LAMP=oC5*<CCRD(9)^(C0R'J(5)^C0«»D(6)/6U. )/60.) 

DO 10 NST=1,5 

1C T(NST)=DIS<STAT(NST,1),STAT(NST,2)) 

DO 11 <sl,9 

11 TQ(.0=eLTA(P)»DELTA(K),T(t<*l)-T(l) 

PRIM 7, TO 

7 FOHmaK'OuS EAST RR s 1930 */» W s •,P9.3,5X,»X s 
:,F9.3,5X,»Z = *,F9.«) 

GOTO 16 
15 STOP 
END 

real function CIS(PFI1,LAH1) 

IMPLICIT REAL(L) 

COMMON PHI, lam 


TIM00050 
TIM00060 
T IM00070 
TIM000«0 
T IM00190 
TO TERMITIMOOIOO 
TIMOOl 10 
TIM00120 
T IM00130 
TIM00190 
TIM00150 
T IM00160 
TIM00170 
TIM00180 
TIM00190 
T 1-00200 

•,F9.3,5x,*Y = •TIM00210 

T IM00220 
TIM00230 
TIMC0290 
TIM00250 
T IM00260 
TM00270 
TIMOO20O 


DATA A/2. 122 2 3 39 35E9/,P/ 2. 12 1099 39 3E9 / , C/9. 966976 7F-1/,F/ 3. 352 3298 TIHO 0290 


69E-3/ 

RETAsAT A\(C*TAN<phI ) ) 

BETA1=ATAU(C*TAN<PpI1)> 

X = ARC0S(SIN(BETA )*SINOETAl )4C0S;BETA) *C0S<5ETAl)*C0S<LAy• 
0=A*X 

P=(A-9)*(X-SIN(X))/(9.*(l.-»C0St''))) 

C=(A-B)*(X4SIN(X))/(9.*(1.-C0S(X))) 

A0=(SIN(BETA)4SIN(PETA1))**2 

RO=(SI\(dETA>-£IN{0ETAl)»**2 

DI «' = ABS (D-AO*P-PO *G ) 

ORIGINAL PAGE IS 

OF POOR QUALITY 


LAMl ) ) 


END 


TIH003C0 
TIM00310 
T 1-00320 
TIM00330 
TIM00390 
TIH00350 
n-00360 
T IM00370 
T IM00380 
TIM0Q390 
TIM00900 
T IM00910 


Figure 3. Time Difference Program. 
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Main Pror.ram 

BP.TA ( 1 ) 

BETA ( 2 ) 

BETA ( 3 ) 

BETA (4) 
CORDd ,2,3) 
C0RD(4,5,6) 
DELTA ( 1 ) 
PELTA(2) 
PELTA(3) 
DELTA (i») 

LAMR 

PHIR 

RPB 

STATd ,1 ) 
:?TAT(2,1 ) 
mT(3,1) 
BTAT(i»,1 ) 
r5TAT(5,1 ) 
1TAT(1 ,2) 
STAT( :,2) 
•=?TAT(3,2) 
BTAT(4,i;) 
STAT(5,2) 

T 

TP 


Baseline - Carolina Beach to Jupiter (W) 
Baseline - Carolina Beach to Cape Race (X) 

Baseline - Carolina Beach to Nantucket (Y) 

Baseline • Carolina Beach to Dana (Z) 

Latitude of receiver; de^ree-mlnute-second 
Longitude of receiver; der;ree-niinute-second 


Secondary Coding Delay at Jupiter (W) 
Secondary Coding Delay at Cape Race (X) 
Secondary Coding Delay at Nantucket (T) 
Secondary Coding Delay at Dana (Z) 


Decimal value of longitude 
Decimal value of latitude 
^IBO 

Latitude of Master (Cape Fear) 
Latitude of Secondary W (Jupiter) 

Latitude of Secondary X (Cape Race) 

Latitude of Secondary Y (Nantucket) 

latitude of "Secondary Z (Dana) 

Longitude of Master (Cape Fear) 
Longitude of Becondary W (Cane Race) 

Longitude of Becodary X (Jupiter) 

Longitude of Secondary Y (Nantucket) 

Longitude of Secondary Z (Dana) 

Distance from receiver to each of the stations 

Time differences for all four pairs 

Figure 4. List of Variables. 
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3ubpro>>ram 


A 

AO 

B 

BETA 

BETA1 

BO 

C 

D 

DIB 

F 

LAM 

LAM1 

P 

PHI 

PHI1 

Q 

X 


Beml-major radius of earth (microseconds) 
Used in curvature correction 
Bemi-minor radius of earth (microseconds) 
Angle of projected point on sphere 
Angle of projected point on sphere 
Used in curvature correction 
B / A 

Spherical arc length 
Geodesic arc length 
Flattening of the earth, 1/29B.2 
Longitude of receiver 
Longitude of Loran station 
Used in curvature correction 
Latitude of receiver 
Latitude of Loran station 
Used in curvature correction 
Angr.lar arc length 


Figure 4 (Cent.). List of Variables. 
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utrd In thii program. The program operatri ot followt: The known valurt for all the 
itation coordinate! of the U.S. Eo«t Coast chain ore initialized at the beginning of the 
program along with the known values of the baseline lengths and coding delays. All 
distances in this program are in microseconds and angles are in radians. When the program 
runs, it requests the entry of the geographies coordinates in the standard degree-mi nute- 
second format. These coordirsotes are then converted to decimal numbers in radians. 
Subprogram "DIS" is called five times to calculate the ore length from the receiver to the 
master and the four secondary stations. These values are then used to solve the time 
difference equation and the time differences just computed are printed out for the four 
pairs. Although this program is set up for the U.S. Eost Coast chain (9930), it may be 
changed to other repetition rates by changing the coordinates of the stations, the base- 
line lengths, and the coding delays. Information on the earth used in this program is 
derived from the Fischer Spheroid System 


IV ADDITIONAL CONSIDERATIONS 

To get a more accurate prediction of Loran-C time differences at a location, there 
are two factors which should be taken into account. These arc: a more accurate model 
of the earth's curvature, and the effect of signal propagation delays. For absolute work, 
there are methods of making estimates of the eorth's shape using techniques such as 
oblate spheroids with unequal axes or power series expansions of the surface. These methods 
ore usually of academic Interest; they are too complicated for a small computer system in 
an aircraft. Using the techniques discussed earlier allows time difference colculctinns to 
be made to within three to four microseconds. 

Propagation delays are difficult to model in o general sense. The terrain over which 
the signals must travel have a slgnificorit effect on lelays.^®^ Many time difference 
calculations ore done assuming a vacuum ond then corrected for salt water retardation. 

For propagation over land, a profile of the surface features may be made, and a typical 
retardation value for the profile may be calculated and applied to subsequent time difference 
calculations . By doing this, accuracy may be improved by two to three microseconds over 
a vacuum depending on the distances from the transmitters and the surface profile used. 

Other factors which can either retard or accelerate signals are atmospheric and ionospheric 
anomalies and local weather conditions. These would have to be monitored daily if absolute 
accuracy is to be obtained. 

Many studies have been made of the effects of the earth's curvoture and propagation 
delays on the accuracy of predicting Loran time differences. Various levels of corrections 
may be applie d depending on the individual sltuatlort. 


V CONCLUSIONS 

The foregoing has presented a technique which is simple enough fo be used on small 
computer systems, yet provides sufficient occurocy for most Loran-C work. \ voriaMms of 
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thr compotrr pro()rom prrirntev^ ha% I'prii Implrmrnfpd tucceufi'lly cw a »moll honv^-hcll^ 
calculator f^Klf>>oypH utinp rlrrnrntary mcKtrU for cursofurp ond propopatioo dplay«, 
tirprp Is pnoupK fIpKibllity to allovs morp occuratp ccwipctiom to bp modp, Thp propram 
discusspd has bppn a prpot aid to rpcpivpr pprformoncp worl^ arsd thp tpstirsp of posltlon- 
flKlrsp calculations . 
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